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TECHNICAL PAPER 


EXACT STATE RECONSTRUCTION IN DETERMINISTIC 
DIGITAL CONTROL SYSTEMS 

I. INTRODUCTION 


Books on modem digital control systems usually address the problem of controlling a continuous- 
time plant driven by a zero-order-hold with a sampled output as shown in Figure 1 [e.g. , Ref. 1 , p. 126]. A 
common solution to this problem is to reconstruct the state of the system at the sampling instant using a state 
observer and then to feed back the reconstructed state [1, p. 195]. However, the state observer has two 
undesirable characteristics. First, it is a dynamical system in itself and, hence, adds additional states and 
eigenvalues to the system which can affect system stability. Second, as a consequence, the reconstructed 
state is normally an approximation to the true state and is usually not a good one early in the state recon- 
struction process. This paper presents a new state reconstructor for deterministic digital control systems 
which has neither of these problems. The new state reconstructor adds no new states, eigenvalues, or 
dynamics to the system. In fact, it does not affect the plant equation for the system in any way; it affects 
only the measurement equation. Furthermore, if the plant parameters are known exactly, the output of this 
new state reconstructor exactly equals the true state of the system. For these reasons, this new state recon- 
structor is herein called the ideal state reconstructor. Useful in the development of this ideal state recon- 
structor are some results to date for continuous-time plants driven by a zero-order-hold with sampled 
outputs. These are reviewed in Section II, prior to the development of the ideal state reconstructor presen- 
ted in Section III. Section IV describes two methods for choosing the reconstructor parameters. Section V 
presents an example illustrating the procedures to completely design the reconstructor using the methods 
described in Section IV. Section VI contains the conclusions and final comments. 
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Figure 1 . Continuous-time plant driven by a zero-order-hold 
with instantaneous measurements. 





II. PRELIMINARY 


For the plant in Figure 1 , x(t) is an n x 1 state vector, u(k) is an r x 1 control input vector, yi(k) 
is an m x 1 output or measurement vector, F is an n x n system matrix, G is an n x r control matrix, 
and Ci is an m x n output matrix. Since yi(k) = Q x(k), where k is the usual shorthand notation for time 
kT, y T (k) represents an instantaneous measure of the system at the sampling instant kT. Hence, the plant 
in Figure 1 can be regarded as having instantaneous measurements for outputs. It is well known that this 
system can be modeled at the sampling instants by the discrete state equations [1, p. 126] 


x(k+l) = Ax(k) + Bu(k) (1) 

y T (k) = Q x(k) , (2) 

where 

4>(t) = [(sI-F)' 1 ] , (3) 

A = <f>(T) , (4) 

and 

T 

B = [ J <j>(\) dx] G . (5) 

0 


4>(t) is the n x n state transition matrix. A and B are the n X n system matrix and then x r control 
matrix, respectively, for the discrete state equations (1) and (2). 

A and B can be determined analytically using equations (3) through (5). An alternative approach, 
which is also quite suitable for numerical computation, is as follows [2]: <|>(t) and /' 4>(\)d\ can be ex- 
pressed in the form of matrix exponential series as 0 


00 



( 6 ) 


and 


2 



( 7 ) 


/V> ax = i ££ . 

0 i = o u u ‘ 

respectively. From equations (6) and (7), 

<t>(t) = I + F [ J 4*X)dx] , 

0 


(8) 


where I is an n x n identity matrix. Hence, J ()>(X)dX can be determined using equation (7) with t = T 
and this result substituted into equation (8) to get 4>(T). With these results, A and B can be found using 
equations (4) and (5). 

Now consider the plant in Figure 2, which is a generalization of the one in Figure 1 . In addition to 
the instantaneous measurement vector yi(kT), the plant in Figure 2 has the measurement vector y F '(kT) 
generated as follows. First, the continuous-time output z(t) is sampled every T/N seconds. Every N 
samples are multiplied by the weighting matrices Hj, j = 0,1,..., N-1 , and then summed to generate the 
output yp(kT) every T seconds. Functionally, this is equivalent to passing the discrete measurements 
generated every T/N seconds through a multi-input/multi-output moving average (MA) process with 
coefficient matrices Hj, j = 0,1 ,...,N-1 [3]. The output of the MA prefilter is sampled every T seconds 
to generate y F (kT). Then y F (kT) has subtracted from it E_ u[(k-l)T], where E_ is a constant matrix, to 
produce the modified MA-prefiltered measurement vector y F '(kT). In Figure 2, C F is a p x n output 
matrix and z(t) is a p x 1 vector. The weighting matrices Hj, j = 0,1,..., N-1 are each q x p. Hence, 
y F (kT) and y F '(kT) are q x 1 vectors. Since u[(k-l)T] is an r x 1 delayed input vector, E_ is a q x r 
matrix. 



Figure 2. Continuous-time plant driven by a zero-order-hold with instantaneous 
and modified MA-prefiltered measurements. 


3 



From Figure 2, y F (kT) can be written as 


z(kT) 

y F (kT) = [H] z(kT-X) 

zfkT-(N-l) X] 

where H is a q x (Np) matrix given by 
H = [H 0 ; H, ; ... ! H n _,] . 

Previously, Polites [4] showed that when 
E_ = HP , 

where fj is the (Np) x r matrix 


C F [J* 4>(k)d\] G 
0 

-(T/N) 

P = C F [J 4>(X)d\] G 
o 

: f -(N-l)(T/N) 

C F LJ cf>(X)d\J G 

L 0 


the discrete state equations for the plant in Figure 2 become 
x(k+l) = Ax(k) + Bu(k) 

x(k) 



I 


I 

I 

j where D_ is a q x n matrix given by 


; D_ = Ha 


( 15 ) 


and a is the (Np) x n matrix 



■ Cp <M0) 

r 

^ a = 

c F <«- i) 


.C,*[-(N-l)iJ . 


(16) 


E_ and D_ can be evaluated analytically using equations (3), (10) through (12), (15), and (16). An 
alternative approach, which can be either analytical or numerical, is as follows. Let t = -j(T/N), where 
j = 0,1 ,.. .,N-1 , and use equation (7) to determine /^ T/N ' <{)(X)dX, j = 0,1 ,...,N-1 . Substitute these results 
into equation (8) to get 4>[-j(T/N)], j = 0,1,...,N-1. At this point, E_ and D_ can be found using 
equations (10) through (12), (15), and (16). 

III. THE IDEAL STATE RECONSTRUCTOR 


A general block diagram of the plant and the ideal state reconstructor is shown in Figure 3 . The 
state reconstructor in Figure 3 is a generalization of the one presented in Reference 5. Observe the 
similarities, and differences, between Figures 2 and 3. In Figure 3, if E_ is given by equation (1 1), then 

x(k+ 1) = Ax(k) + Bu(k) (17) 

and 


y F '(k) = D_ x(k) , (18) 

where D_ is given by equation (15). This follows from Figures 2 and 3 and equations (13) and (14). Also, 
in Figure 3, 


y F "(k) = (DJ D_)-> D_ t y F '(k) 


(19) 
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PLANT 


IDEAL STATE RECONSTRUCTOR 



Figure 3. General block diagram of the plant and the 
ideal state reconstructor. 


However, for equation (19) to be meaningful, (D_ T D_) ' must exist, and this occurs only when (D_ T D_) 
is nonsingular. Recall that D_ is a q X n matrix. If q 2= n and D_ has maximal rank (i.e., rank n), then 
(D_ T D_) is positive definite and therefore nonsingular [6]. Hence, equation (19) requires that q S 5 n and 
rank (D_) = n for it and the ideal state reconstructor to be meaningful. Assuming this is the case, it 
follows from equations (17) through (19) that the discrete state equations for the system in Figure 3 are 


x(k+l) = Ax(k) + Bu(k) 


( 20 ) 


y F "(k) = x(k) . 


Hence, the output of the state reconstructor, y F "(k), equals the true state, x(k), exactly. Moreover, this is 
achieved without adding any new states, eigenvalues, or dynamics to the plant, since the plant equation 
(20) for the system in Figure 3 is identical to the plant equation ( 1 ) for the plant in Figure 1 . 

In summary, the requirements which must be satisfied for the ideal state reconstructor in Figure 3 
to exactly reconstruct the state of the plant in Figure 3 are as follows. The plant matrices F, G, and C F 
must be known exactly in order to determine E_ and D_ according to equations (11) and (15), respec- 
tively. The number of rows, q, in the weighting matrices Hj, j = 0,1,...,N-1 must be greater than or 
equal to the number of states in the plant, n. The qxn matrix D_ in equation (15) must have rank n. So 
long as these requirements are met, there are no other restrictions on the ideal state reconstructor, includ- 
ing its weighting matrices. This is very desirable because it means there are a multitude of ways to choose 
the weighting matrices and still achieve exact state reconstruction. Two distinct methods for doing this 
are presented in Section IV. 
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IV. CHOOSING THE PARAMETERS IN THE IDEAL 
STATE RECONSTRUCTOR 


Once the measurements which are inputs into the ideal state reconstructor have been specified, 
the p x n output matrix C F is defined. To complete the definition of the reconstructor, the following 
parameters must be chosen: the number of weighting matrices in the state reconstructor N, the number of 
rows in each weighting matrix q, and the elements of each q x p weighting matrix Hj, j = 0,1,...,N-1. 
While there are countless ways of choosing these, two distinct methods are presented. 


Method A. Let H = (o^a)' 1 a T 

Consider H, the q x (Np) matrix formed by catenating the weighting matrices Hj, j = 0,1,..., N- 
1 as in equation (10). Consider also a, the (Np) x n matrix defined by equation (16). If Np 5 n and 
rank(a) = n, then (a T a) is positive definite and therefore nonsingular [6], If this is the case, then H can 
be given by the pseudo inverse of a or 


H = (a T a)- ] a T 


( 21 ) 


Hence, N must be chosen so that Np ^ n, or equivalently N 5= n/p. Having chosen N, a can be deter- 
mined using equation (16). If it has rank n, then let H be given by equation (21). This makes H an 
n x (Np) matrix and means that q = n in this case. Substituting equation (21) into equation (15) yields 
D_ = I, where I is an n x n identity matrix. Consequently, 


(D_ T D_)"' D_ t = I 


( 22 ) 


also. This simplifies the state reconstructor by eliminating the need to solve equation (19) in it. From 
equations (11) and (21), E_ is found to be 


E_ = (a T a)'' a T (3 , 


(23) 


where (3 is defined by equation (12). From equations (10) and (21), 


[H 0 H, 


H n _,] = (a T a)-' a T 


(24) 


which reveals the N n x p weighting matrices. The state reconstructor is now completely defined for 
this method. It is equivalent to the one described in Reference 5. 
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Method B. Let H = I 


In this method, let 


H = I 


( 25 ) 


where I is an (Np) x (Np) identity matrix. Since H is a q x (Np) matrix in general, this means that 
q = Np in this case. Recall from Section III that one requirement for the reconstructor to be meaningful 
is for q ^ n. Consequently, N must be chosen so that Np 5= n or, equivalently, N ^ n/p, just as in 
Method A. Having chosen N, a can be determined using equation (16). Substituting equation (25) into 
equation (15) yields 


D_ = a 


(26) 


Recall also from Section III that another requirement for the reconstructor to be meaningful is for 
rank(D_) = n. By virtue of equation (26), it is necessary for rank(a) = n, just as in Method A. It 
follows from equation (26) that 


(D_ T D_)-' D_ t = (a T a)-' a T . 


(27) 


From equations (21), (22), (25), and (27), one can see the relationships which exist between H and 
(D_ t D_) _1 D_ t in Methods A and B. From equations (11) and (25), E_ is found to be 


E- = (3 


(28) 


where (3 is defined by equation (12). From equations (10) and (25), 


[Ho ! H, 


H n .,] = I 


(29) 


which reveals the N (Np) x p weighting matrices. However, recognize from equations (9) and (25) that 


yp(kT) = 


r ?(kT) 


?( kT - f ) 


L z[kT - (N-l) -2j-) 


(30) 
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Hence, in this method, the intermediate output y F (kT) can be formed directly by catenating the meas- 
urements z[kT - j(T/N)], j = 0,1,..., N-l as in equation (30). This is the beauty of this method. 

This description completes the definition of the ideal state reconstructor for this method. In Sec- 
tion V, an example is presented which illustrates the procedures to completely design the ideal state 
reconstructor using both methods presented in this section. 


V. AN EXAMPLE 


Consider the double integrator plant driven by a zero-order-hold as shown in Figure 4. The con- 
tinuous-time output z(t) is sampled every T/N seconds and input into the ideal state reconstructor together 
with the control input u(kT). The parameters in the state reconstructor will be chosen to achieve exact 
state reconstruction using both of the methods described in Section IV. 

Manipulating the plant in Figure 4 into the format of Figure 3 yields 


[F] 


0 1 
0 0 


(31) 


[G] 



(32) 


and 


[C F ] = [1 0] . 


(33) 


Using equations (31) and (32) and the formulas presented in Section II, 


4>(t) = 


1 t 

0 1 


f <KX)d\ = 
o 


t 


0 


(34) 


(35) 
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Figure 4. Plant and ideal state reconstructor 
for the example. 


A 



T 

1 


and 


B 



Since Cp is a p x n matrix, it follows from equation (33) that p = 1 and n = 2. Hence, the requirement 
to select N so that N s* n/p can be satisfied by letting 


N = 4 . 

Using equations (12), (16), (32) through (35) and (36), a and (3 are found to be 


(36) 


a 


1 

1 

1 

1 


0 

T 

T 

2T 

T 

3T 

T 


(37) 
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and 


0 


T 2 



4T 2 

~32 

9T 2 

IT 


(38) 


respectively. In equation (37), eliminating any 2 rows forms a 2 x 2 matrix with nonzero determinant, 
assuming of course that T > 0. Hence, rank(a) = 2 = n and so (a T a) is nonsingular. Consequently, 
(cJa)" 1 a T exists and is found to be 


(ot T a)-' a T 



(io) (io) 

(w) (4) 



(39) 


using equation (37). 


Method A. Let H = (a T a) _1 a T 

From equations (21), (24), and (39), 






H = [H 0 H, 1 H 2 ! H 3 ] = 





which reveals the weighting matrices Hj, j = 0,1, 2, 3. From equation (22) and the fact that n = 2, 


(D_ t D_)-' D_ t 


1 0 
0 1 
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From equations (23), (38), and (39), 


E_ = 



This completely defines the ideal state reconstructor for this example and this method. 

Method B. Let H = I 

Since N = 4 and p = 1, it follows that q = Np = 4. From equations (25) and (29), 


H = [H 0 1 H, | H 2 ; H 3 ] = 


0 

0 


1 

0 


0 

1 


0 

0 


which reveals the weighting matrices Hj, j = 0,1, 2, 3. However, from equation (30), the intermediate 
output yp(kT) can be related directly to the measurements z[kT - j(T/4)], j = 0,1, 2, 3 and is 


yp(kT) = 


z(kT) 


Z (kT--) 
z(kT — ~ ) 

4 

z(kT-^) 

4 


From equations (27) and (39), 


(D 1 D_) 1 D_ t = 


'(b) 

(w) 


U) (to) (to) 
(st) (-5t) (^) 


J 
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From equations (28) and (38), 


E_ = 


0 

T 2 

~32 

4T 2 

"32" 


9T 2 

~32 


This completes the definition of the ideal state reconstructor for this example and this method. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


A new state reconstructor for deterministic digital control systems has been presented which 
offers two distinct advantages over the widely used state observer. First, it adds no new states, eigen- 
values, or dynamics to the system and, consequently, will not alter the stability of the system. In fact, 
adding the new state reconstructor to the system will not affect the plant equation in any way. It only 
affects the measurement equation. Second, if the plant parameters are known exactly, the reconstructed 
state will exactly equal the true state of the system, not just approximate it. For these reasons, this new 
state reconstructor has herein been called the ideal state reconstructor. Its disadvantages are that the plant 
output must be sampled and calculations must be performed more frequently than with the state observer. 
Fortunately, this is not the problem it was 20 years ago, considering the speed of today’s digital 
computers. 

If the research in this report is extended, two approaches are recommended. One is to explore 
other methods for choosing the parameters in the ideal state reconstructor. Two distinct methods have 
been presented here, each with attractive features of its own. However, these do not exhaust all the possi- 
bilities, by any means. For example, the weighting matrices could be selected so the MA prefilter acts as 
a multi-input/multi-output low-pass filter for the case where measurement noise is present. The other 
approach is to investigate the robustness of the ideal state reconstructor and see how it compares with the 
state observer’s. Specifically, the following questions should be addressed. What effect do modeling 
errors in the plant have on the ideal state reconstructor, and how does this compare with the state 
observer? What effect do plant process and measurement noise have on the ideal state reconstructor and 
how does this compare with the state observer, or even the Kalman filter? How can the robustness of the 
ideal state reconstructor be improved? Increasing the number of weighting matrices, N, may be one 
possibility. Catenating the ideal state reconstructor with a state observer, or a Kalman filter, may be 
another. This might produce a composite estimator which is better than either the ideal state reconstruc- 
tor, the state observer, or the Kalman filter alone. 


13 


REFERENCES 


1. Jacquot, R. G.: Modern Digital Control Systems. Marcel-Dekker, New York, 1981. 

2. Cook, G.: Lecture notes from EE312 - Digital Control Systems, Vanderbilt University, Fall 1984 
(Personal Notes). 

3. Goodwin, G. C., and Sin, K. S.: Adaptive Filtering, Prediction, and Control. Prentice-Hall, 
Englewood Cliffs, New Jersey, 1984, p. 32. 

4. Polites, M. E.: Modeling Digital Control Systems With MA-Prefiltered Measurements. NASA 
TP-2732, George C. Marshall Space Flight Center, Huntsville, Alabama, June 1987. 

5. Polites, M. E.: A New Approach to State Estimation in Deterministic Digital Control Systems. 
NASA TP-2745, George C. Marshall Space Flight Center, Huntsville, Alabama, July 1987. 

6. Greville, T. N.: The Pseudoinverse of a Rectangular or Singular Matrix and Its Application to the 
Solution of Systems of Linear Equations. SIAM Review, Vol. 1, No. 1, 1959, p. 38. 


14 


2. government accession no. 


1. REPORT NO. 

NASA TP-2757 


4. TITLE AND SUBTITLE 

Exact State Reconstruction in Deterministic Digital 
Control Systems 


7. AUTHOR(S) 

Michael E. Polites 


9. PERFORMING ORGANIZATION NAME ANO ADORESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12. SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

August 1987 


6. PERFORMING ORGANIZATION CODE 


8. PERFORMING ORGANIZATION REPORT U 


10. WORK UNIT NO. 

M-565 


1 1. CONTRACT OR GRANT NO. 


13. TYPE OF REPORT & PERIOD COVERED 

Technical Paper 


14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 

Prepared by Structures and Dynamics Laboratory, Science and Engineering Directorate. 


16. ABSTRACT 

This report presents a new state reconstructor for deterministic digital control systems which is ideal 
in the following sense: if the plant parameters are known exactly, the output of the new state reconstructor 
will exactly equal the true state of the plant, not just approximate it. Furthermore, this ideal state recon- 
structor adds no additional states or eigenvalues to the system. Nor does it affect the plant equation for the 
system in any way; it affects only the measurement equation. While there are countless ways of choosing 
the ideal state reconstructor parameters, two distinct methods are described here. An example is presented 
which illustrates the procedures to completely design the ideal state reconstructor using both methods. 


17. KEY WORDS 

Digital Control Systems 
State Estimation 
State Reconstruction 
State Reconstructor 


18. DISTRIBUTION STATEMENT 

Unclassified/Unlimited 

Subject Category: 31 


19. SECURITY CLASSIF. (of thU roportl 

Unclassified 


MSFC - Form 3292 (May 1969) 


20. SECURITY CLASSIF. (of thU page) 

21. NO. OF PAGES 

Unclassified 

18 



For tale by National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1987 














